Controlled transport of matter waves in two-dimensional optical lattices 
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We propose a method for achieving dynamically controllable transport of highly mobile matter- 
wave solitons in a driven two-dimensional optical lattice. Our numerical analysis based on the 
mean-field model and the theory based on the time-averaging approach, demonstrate that a fast 
time-periodic rocking of the two-dimensional optical lattice enables efficient stabilization and ma- 
nipulation of spatially localized matter wavepackets via induced reconfigurable mobility channels. 
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Controlled manipulation of stable, spatially localized 
matter waves, similar to routing of optical pulses in pho- 
tonic devices, is a very attractive goal from the viewpoint 
of emerging integrated technologies based on the use of 
ultracold atomic gases - Bose-Einstein condensates 
In recent years optical lattices have been suggested as 
an instrument for such a control. Long-distance trans- 
port of ultracold atoms in a one-dimensional Bessel opti- 
cal lattice potential was demonstrated experimentally Q ■ 
Theoretical studies of matter- wave solitons, loaded into 
a rapidly driven asymmetric one-dimensional optical lat- 
tice, have demonstrated that such an "optical ratchet" 
could enable their mass-dependent transport [§]. Trans- 
port of matter waves in a two- or three-dimensional trap- 
ping geometry is a more complex and challenging task, 
especially considering the intrinsic instability of the self- 
localizing condensate with the negative scattering length. 

Nonlinearly localized atomic wavepackets - matter- 
wave solitons, created in a two- or three-dimensional 
trapping geometry, are predicted to be stabilized by 
static optical lattices ■ Other methods of stabilization 
rely on the time-periodic management of the conden- 
sate properties or non-local interaction between atoms 
5} . Stability of the matter waves in two-dimensional op- 
tical lattices and under different conditions of the time- 
periodic management (see, e.g., @) has been extensively 
studied. The main challenge still remains: To suggest 
an efficient method for non-destructive, dynamically con- 
trolled transport of the stabilized atomic wavepackets. 

The problem of transport of localized nonlinear excita- 
tions is acutely posed in several cross-disciplinary areas of 
physics (see, e.g., Ref. Q and references therein) because 
periodic systems, in general, greatly inhibit the mobility 
of the localized states. In particular, mobility of discrete 
solitons is an important problem in nonlinear photonics 
and lattice dynamics [8|, |9|. In periodic structures and 
waveguide arrays, the motion of spatially localized opti- 
cal beams is affected substantially by the lattice period- 
icity, so that the effective Peierls-Nabarro potential |8[ 
induced by the lattice prevents their free motion. It was 
shown that moving solitary waves in discrete lattices re- 
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move only at some limited velocities 
discrete optical solitons in two-dimensional photonic lat- 
tices was also found for the case of quadratic nonlinear 
response that leads to parametric coupling of the funda- 
mental and second-harmonic fields and results in their 
mobility after an initial kick [l2j |. 

In this Letter we propose and analyze the mechanism 
for controlled transport of two-dimensional matter-wave 
solitons, created in a Bose-Einstein condensate of atoms 
with a negative scattering length. The transport is real- 
ized by means of a rocking two-dimensional optical lat- 
tice [13j] , where the term " rocking" refers to time-periodic 
shaking of the lattice 14|. Our analysis, based on the 
mean-field model, demonstrates that the fast rocking of 
the lattice enables dynamical creation of reconfigurable 
mobility channels for matter waves. Similar channels 
were previously studied in the context of lattices that do 
not coincide in dimensionality with the wavepacket and 
guide a matter wave along the fixed free spatial direction 
hj. Here we show that the rocking optical lattices enable 
both the efficient stabilization and dynamical "routing" 
of nonlinearly localized matter wavepackets. 

To describe the dynamics of solitons in a Bose-Einstein 
condensate loaded into a two-dimensional lattice, we con- 
sider a mean-field model of an ultracold atomic cloud in a 
pancake trapping geometry with an optical lattice poten- 
tial aligned with the weak trapping directions. As long 
as the transverse collective modes of the condensate in 
the direction perpendicular to the optical lattice are not 
excited, the system can be treated as two-dimensional 
and described by the Gross-Pitaevskii equation: 

^ + Vi^+HV + ^OL(r,^ = 0, (1) 

where V\ = d 2 /dr 2 , and r = (x,y). This model is 
derived by assuming the units of energy, length, and 
frequency: El = fi. 2 fc|/(2rn), ar, = l/k^, and lot, = 
Er,/h, respectively, where m is the atomic mass, and 
fc_L is the wavevector of the optical lattice. The three- 
dimensional mean-field model is reduced to the two- 
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dimensional equation by assuming that the wavefunction 
is separable: ipa D (r,z) = ip2D(r)<f>iD(z), where 4>id{z) 
is the normalized ground state wavefunction of a one- 
dimensional harmonic potential with the trapping fre- 
quency uj z . With these assumptions, the wavefunction ip 
in Eq. JT]) relates to tp2D as follows: ip = ip2D\/92D, where 
92D = 4v / 7r(a s /ao)(o; z /a;L) 1 ^ 2 is the renormalized coef- 
ficient that characterizes interaction of the condensate 
atoms with the s-wave scattering length a s . The number 
of atoms is given by: TV = N/g 2 D, where N = J \^\ 2 dx 
is the norm of the dimensionless wavefunction. 

For the 7 Li atoms with a s = —0.21 nm, trapped using 
a CO2 laser with A = 10.62 /im [lij], the physical values of 
the characteristic scales are: ao = 845 nm, loq = 2ttx 2039 
Hz. Consequently, the rescaled interaction coefficient is 
cj2D = 1-03 x 10~ 3 , and the typical soliton has the norm 
N s» 10, which corresponds AT ~ 10 4 . 




FIG. 1: (Color online) Norm of the wavehmction correspond- 
ing to a stationary soliton state in a rocking lattice potential 
VoL(r, t) at t = (curve a — b) and in a time-averaged po- 
tential F S (R S ) (curve c - d) with Y = 0, X = 2.4 . (a-d) 
Spatial density structure of the matter-wave solitons at the 
marked points on the top panel. 



The optical lattice potential in Eq. (|T|) is created by 
two pairs of counter-propagating laser beams and has the 
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FIG. 2: (Color online) Dynamics of a moving soliton in the 
rocking lattice obtained by solving Eq. fT]) numerically with 
periodic boundary conditions. Initial states are: (a) a station- 
ary soliton of the potential Voh(r, 0) at /i = —0.7 [Fig. Q] (b)] 
and (b) a stationary soliton of the time-averaged potential 
K(Rs) at Y = 0, X = 2.4, ft = -0.7 [Fig. [^d)]. 



following form: 

VbL(r > t) = V r B cos[a; + A:(t)] + ^ / cos[y + K(t)]. (2) 

It is driven periodically, i.e. the phase detuning be- 
tween the laser beams forming the lattice is modulated 
as: X(t) = X sin(ujt), and Y(t) = Y sin(u)t). By using 
Fourier decomposition, one can see that the optical lat- 
tice has both a static "backbone" and a time-dependent 
component. This distinguishes transport in the poten- 
tial ([2]) from the ratchet-type transport of matter-wave 
solitons 3]: Here the time-average force acting upon a 
wave-packet at any given point in space is non-zero. 

At uj = (or t = 0), the static optical lattice VoL(r, 0) 
supports stationary matter-wave solitons [3] in the form: 
ip(r,t) = ^(r) exp(-fit) at the values of chemical poten- 
tial, /«, below the lower edge of the first spectral band 
(see Fig. [TJ. These are found numerically by solving the 
stationary version of the model Eq. (UJ [17|. According 
to the Vakhitov-Kolokolov stability criterion, such local- 
ized states are dynamically stable away from the band 
edge, where d/i/dN < The mobility of the solitons 
is suppressed due to the energy difference between the 
nonlinear localized states with the same TV but different 
symmetry, that needs to be overcome in order to move 
the initially stationary state across the lattice. This is 
analogous to the Piers-Nabarro potential in discrete sys- 
tems [8j. On the other hand, it is the coupling to the 
lattice that both stabilizes the solitons against the col- 
lapse and allows us to manipulate the localized states 
by changing parameters of the lattice. The typical norm 
of a two-dimensional soliton considered in our dynamical 
simulations below is near the collapse threshold for a free 
soliton, which for our normalization is: N* — 11.7 0]. 

To find the regime of enhanced mobility in the lattice 
we assume that a mobile soliton does not significantly 
change its shape and can be described solely by the dy- 
namics of the center of mass, r (t) = (xo(t),y (t)). This 
assumption is valid provided that the driving frequency, 
uj, is much greater than the eigenfrequency of soliton 
width oscillations (ujq ~ 1 in our system). Following the 
treatment in [19j | , the Hamiltonian theory that treats the 
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soliton as a classical particle in an effective potential [20( 
can then be invoked to derive the equations of motion for 
the center of mass. The matter-wave soliton loaded into 
the lattice is approximated by the Gaussian function: 



# (r(t),r (t)) = Aexp 



(x - x ) 2 + (y - Vof 



2a 2 



The variational theory (see, e.g., [6j) shows that this ap- 
proximation is in good agreement with the exact station- 
ary solutions for values of fi far from the band edge, where 
the soliton is well localized. In this case the soliton's am- 
plitude and width are related as follows: N = irA 2 a 2 , and 
N = 2tt[2 - V a 4 exp(-a 2 /4)]. The effective potential for 
the soliton-as-particle motion is then found as follows: 



1 

V eS {r ,t) = — 



\V (r,r o )\ 2 VoL(r,t)dr, 



where dr = dxdy, and the equations of motion for the 
center of mass, d 2 r /dt 2 = — dV e g (r )/dr , read: 

d 2 x 4 d 2 y _<£ . , v , 

—— = V x e * sm[xo - X ), — j- = V y e * sm(y - Y). 

The motion of the soliton's center of mass can be sep- 
arated into the slow and fast components: Xo(t) — 
X„(t) + £(*), y (t) = Y s (t) + r)(t), where the typical 
evolution of the slow variables R s = (X S ,Y S ) occurs on 
the time scale r 3> a; -1 . The time averaging procedure 
Til 21 1 is then applied to derive equations of motion for 
the slow (compared to uS) dynamics of a soliton center 
of mass in a driven lattice: d 2 R s /dt 2 = — <fl4(R s )/<iR s . 
The time-averaged effective potential is separable: 



V S (R S ) = A sx cos(2X s ) -(- A sy cos(2r 5 .), 



(3) 



with the strength explicitly depending on the driving fre- 
quency, u>, and on the characteristic width of the soliton 
loaded into the lattice, a: 



V 2 

A — x r~ 

sx ~ So; 2 ' 
V 2 



{4[J 1 (Xo)f-[J 2 (X )] 2 }, 

A sy = ^e-^{4[j 1 (y )] 2 -[j 2 (y )] 2 }, 
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where J n are Bessel functions of the first kind. We stress 
that the time-averaged potential is conservative even in 
the presence of driving at certain values of (Xq,Yq), 
namely such that Jo(X ) = or Jo(Y ) = 0. 

The time-averaged conservative potential, V s , supports 
stationary matter-wave soliton solutions depicted by the 
family (c-d) in Fig. [TJ These solitons may be dynamically 
stable in the entire existence region of the solitons sup- 
ported by the static lattice potential VoL(r, 0) [HI due 
to the fact that the band edge in V s is shifted to lower 
values of the chemical potential, [i (dashed line Fig. [TJ. 
This indicates that the driving may stabilize a moving 



solitons below the collapse threshold, as predicted in [6( 
for amplitude-modulated optical lattices. In addition, 
for the same values of /j,, the mobility of the solitons sup- 
ported by V s is expected to be enhanced compared to 
those in the static lattice VoL(r, 0) due to the selectively 
suppressed spatial modulation. 

In order to confirm the predictions of improved mobil- 
ity of the localized states in the driven lattice, we com- 
pare the dynamics of two initially stationary states cor- 
responding to /j, = —0.7 and supported by Vol(i",0) and 
V s potentials (see Fig. At the time t = both the 
initial states are given exactly the same initial velocity 
in the x-direction. Nevertheless, the soliton, supported 
by the static optical lattice Vol(i",0) (point b in Fig. [TJ) 
decays almost instantly, whereas the soliton supported 
by the V s in the driven lattice (point d in Fig. [TJ) moves 
and retains its shape [Fig. [2] (b)]. 
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FIG. 3: (Color online) Dynamical control of the optical lat- 
tice parameters enabling creation and switching of the soliton 
mobility channels in the time-average potential. Left column: 
Rocking optical lattice potential, VoL(r, i) at the instances of 
time t\ and ti\ the lines marking x = and y = are a guide 
to the eye. Middle column: Time dependence of the rocking 
amplitudes. Right column: Time-averaged lattice potential, 
V s (R.s), corresponding to Yo = 0, A sy = 0.5, A sx — 0.002 (at 
h) and X a = 0, A sx = 0.002, A sy = 0.5 (at t 2 ). 

Next, we explain the basic principles of dynamical 
transport control of the mobile solitons in the driven lat- 
tice. While the lattice potential, Vol(i", f) is not static, 
and the positions of its maxima and minima are differ- 
ent at various times (see Fig. left column), its action 
on the soliton center of mass during the periods of time 
t oj^ 1 can be described by the conservative potential 
V s . By examining the form of the time-averaged poten- 
tial ©, one can see that the amplitude of the spatial 
modulation can be effectively suppressed in one of the 
orthogonal directions by controlling the amplitude of the 
periodic driving. The most dramatic example of such 
suppression is presented in Fig. [3l During the periods 
of time when the amplitudes of driving Xq and Yq are 
kept constant (segments marked by the points t\ and t% 
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FIG. 4: (Color online) (a,b) Center of mass trajectory of a 
moving soliton corresponding to \i = —2.0 in a rocking lattice 
potential VoL,(r, t), superimposed onto the contour plot of 
the potential at t = 0, and, in (c), onto the density profile 
of the moving soliton. Sharp turning points correspond to 
the switching of the mobility channels as shown in Fig. (d) 
The profiles of the moving soliton corresponding to the points 
1 (solid), 2 (dashed), and 3 (dash-dotted) in (c). 



perturbed by small amplitude fast oscillations. By ex- 
amining the cross-sections of the soliton profile, we also 
observe that the localized state is only weakly affected 
by these dynamical manipulations [see FigHJd)], which 
justifies our soliton-as-particle approach. 

In conclusion, we have demonstrated, both analytically 
and numerically, that time-modulated optical lattices en- 
able realization of clean, controlled transport of matter 
waves in more than one spatial dimensions. In doing 
so, we have used physical parameters relevant to exper- 
iments with the ultracold atoms. This work opens up 
a fascinating possibility to realize transport and rout- 
ing of matter-wavepackets in lattice potentials of more 
complex symmetry, as well as in fully three-dimensional 
trapping geometry. Additional challenge is to achieve 
the transport of gap solitons supported by a lattice in 
a condensate with a positive scattering length, that are 
strongly "pinned" by the lattice and typically immobile. 
The work along these directions is currently underway. 

This work was partially supported by the Australian 
Research Council (ARC). The authors are grateful to Dr. 
T. Alexander for stimulating discussions. 



in Fig. [3J middle column), one of the amplitudes, A or 
Y , is set to a value such that the time-averaged potential 
has a form of a two-dimensional lattice with the spatial 
modulation amplitude A sx or A sy radically suppressed, 
but not equal to zero (see Fig. [3J right column). Thus 
the mobility channels for a localized state are effectively 
created in the two-dimensional lattice along one of the 
orthogonal directions. This effective lattice is fully re- 
configurable in real time, since direction of the mobility 
channels can be switched as shown in Fig. [3J During 
the switching time, which is short compared to u> , the 
concept of time-averaged potential is not applicable. 

To demonstrate the transport of a matter-wave soli- 
ton in the dynamically reconfigurable driven lattice, we 
prepare the initial state in the form of the exact sta- 
tionary soliton supported by the potential V s [point (c) 
in Fig. [Tj. We then load it into the lattice © driven 
with the frequency uj = 10 and the amplitudes X a = 2.4, 
Yq = 0, and let it evolve. At the time t « 60 the driving 
amplitudes are rapidly switched to the values Xq = 0, 
Yq = 2.4, as shown in Fig. [3J (middle panel). The re- 
sulting trajectory of the soliton's center of mass is shown 
in Fig. 2Ja) against the contours of the potential Vol at 
t = 0. The trajectory is easily reconfigurable by changing 
the sequence of switching [see Fig. HJb)]. A more com- 
plex routing of the soliton along orthogonal directions in 
the lattice can also be realized, as shown in Fig. Etc). 
In agreement with our initial assumption, the center-of- 
mass trajectory consists of slow (compared to uS) motion 
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